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We study theoretically the behavior of thermal massless Dirac fermions inside graphene-based Ferromagnetic | 
Insulator | d-wave/s-wave superconductor (F|I|d and F|I|S ) junctions in the ballistic regime. Using the Dirac- 
BdG wave functions within the three regions and appropriate boundary conditions, the Andreev and Normal 
reflection coefficients are derived. By employing the obtained Andreev and Normal reflection coefficients 
the characteristics of heat current through the F|I|d and F|I|S junctions are investigated within the thin 
barrier approximation. We find that for s-wave superconductors, thermal conductance oscillates sinusoidally 
vs barrier strength. The finding persist for the values of a, the orientation of d-wave superconductor crystal 
in the fc-space, below 7r/4. By increasing temperature, the thermal conductance is increased exponentially 
for small values of a and for larger values the quantity is modified to exhibit a linear behavior at a — n/ 4 
which is similar to Wiedemann-Franz law for metals in low temperatures. 



I. INTRODUCTION 

Graphene is a single layer of carbon atoms which was 
introduced to the scientific community by Novoselov et 
al. in 2004^ ^. Most of applicable and interesting char- 
acteristics of graphene has been investigated intensely 
experimental!}' ^'^ . Because of interesting phenomena 
which graphene showed, the artificial material has been 
received robust attentions from theoretical and experi- 
mental physics communitieJ2HIS!_ Induction of supercon- 
ductive correlations into graphene layer by proximity of 
superconducting electrodes observed experimentally by 
Heersche et al^. Also inducing the ferromagnetism into 
graphene layers by means of proximity effects observed 
experimentally by Tombros et al^^. In this regard, 
theoretical scientists utilized Dirac-Bogoliubov-de Genne 
(Dirac-BdG) for investigating and predicti ng inte resting 
phenomena because of the proximity effectJ^I^IinHIll^ 

Interplay between ferromagnetic graphene sheets 
and conventional superconductor is generalized 
theoreticall}^. Also Linder et al^ generalized the 
theoretical investigations from conventional to uncon- 
ventional superconductors in similarity with metallic 
cases which has been studied intenseljff^llSSl . Most of 
the previous works are devoted to study electronic trans- 
port properties of the graphene-based junctions such 
as; Josephson currents^ electronic conductanc^^Ell^ 
spintronic conductanc e^^'^"^' , shot nois d^"*^^ and etc., 
but poor attentions has been focused on the heat 
transport properties and electronic thermal conductance 
of the junctions. The BTK formalism is generalized 
by Bardas and Averirl^^ for obtaining electronic ther- 
mal conductance in the clean limit. Also, Devyatov 
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et alW^ studied electronic thermal conductance of 
Normal metal|I|d junctions in the ballistic regime. 
For high-sensitive devices including graphene-based 
junctions, knowing all electronic and thermal proper- 
ties of the junctions are crucial and important points 
from application point of view. Few previous works 
devoted for investigating of electronic thermal trans port 
characteristics of the graphene-based j unction J^HISI]. 

In this paper we especially investigate signatures of 
(i2;2_j^2-wave symmetry on the electronic heat transport 
characteristics of the F|I|d junctions. We start with the 
Dirac-BdG Hamiltonian and use obtained wave func- 
tions within the three regions and appropriate bound- 
ary conditions at interfaces for deriving the Andreev 
and Normal reflec tion coefficients in the thin barrier 
approximatiorP^m, Using the mentioned coefficients we 
present numerical investigations of electronic thermal 
conductance of the Ferromagnetic|Insulator|(i/s-wave su- 
perconductor junctions in the ballistic limit. We find that 
for F|I|d junctions the electronic thermal conductance, F 
shows oscillatory behavior vs barrier strength. Increas- 
ing the orientation angle of the superconducting gap a 
up to values close to 7r/4 only enhances whole values of 
the F and for maximum value of the superconducting 
gap orientation angle namely a — tt/A magnitude of the 
oscillations diminish highly. Another finding is that the 
electronic thermal conductance shows an exponential in- 
crease vs temperature for small values of a. By increasing 
the angle of superconducting gap orientation up to 7r/4 
the mentioned exponential form modify to linear form 
and at a — n/A: the thermal conductance shows precisely 
linear increase with respect to temperature, namely F cx T 
that the finding induces in mind the Wiedemann-Franz 
law^^ from metals in the low temperatures. The paper is 
organized as follows: 

In Sec. [njwe explain the analytical derivations of An- 
dreev and Normal refiection coefficients by starting from 
the Dirac-BdG Hamiltonian and in Sec. Illll the electronic 
thermal conductance of F|I|S is investigated by plotting 
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the quantity with respect to strengthes of exchange field 
and barrier in the thin barrier approximation regime. In 
Sec. [TVlwe study the effects of d^2_^2-wave symmetry 



on the electronic thermal conductance of the junctions 
within the thin barrier approximation. The paper go to 
end with conclusions in Sec. Ivl 
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II. THEORY 

We study interplay between graphene-based ferromag- 
netic and superconductor junctions in the ballistic limit, 
therefore we employ the Dirac-BdG Hamiltonian for ob- 
taining suitable wave functions. The general Dirac-BdG 
equation incorporating ferromagnetism and supercon- 
ductivity reads as^'^: 



Hq — ah 
A*(T) 



A(r) 

-{Ho -ah) 



Ho{r) = -ihvpicxdx + 'yydy) + U{r) - Ep (1) 

where ax and ay are 2x2 Pauli matrices and A(T) 
stands for temperature-dependent order parameter of su- 
perconducting region, also h represents the strength of 
exchange field in the ferromagnetic region, eg- stands 
for excitation energies of holelike and electronlike quasi- 
particles. For obtaining Dirac-BdG wave functions in 
the Normal, ferromagnetic and superconducting region, 
one should set {h = 0,A(T) = 0}, {h 0, A(T) = 0} 
and {h — 0, A ^ 0} inside Eq. ([!]), respectively. Here 
a — ±1 stands for spin- up and -down quasiparticles and 
a = —a. Also U{r) shows the Fermi energy mismatch. 
Throughout the paper we consider a step function for 
spatial-dependency of the superconducting gap, namely 
A(a;,r) = A{T)e{x) in which e{x) is the weU known 
step function. By solving the Eq. ^ in the Ferromag- 
netic region we obtain Dirac-BdG wave functions for elec- 
tronlike and holelike quasiparticles as follows: 
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where 6a{6'g) are propagation angles of electronlike 
(holelike) quasiparticles with respect to the normal tra- 
jectory into the interface at a; = 0. We define the two 
incident angles as 



if, — arcsm 
'ct = arcsin | 



e+Ep+rrh , 

hvpq 
e—Ep+ah 



(3) 



and x-components of the wave vectors for electronlike 
and holelike quasiparticles in the Ferromagnetic region 
are obtain as 



J, _ t-Ep+ah 



COS y^, 



(4) 



FIG. 1. (Color online) Model for suggested set up 
of FerromagneticjBarrier|ci-wave superconductor graphene- 
based junctions. The junction is located at a; = 0. Axes 
of d-wave superconductor crystal are represented with a and 
b. The superconductor crystal orientation angle with respect 
to normal trajectory (x-axis) is shown by a. 



in which Ep and q are Fermi energy and y-component 
of wave vector in the Ferromagnetic region, respectively. 
For normal region with a barrier potential Vq, it is suf- 
ficient to set h = and Ep [Ep — Vq) in the above 
obtained equations inside the Ferromagnetic region. The 
barrier potential Vq can be applied by a gate voltage 
into the region. Within the d-wave superconductor re- 
gion (a: > 0), Dirac-BdG wave functions for electronlike 
and holelike quasiparticles are obtained: 



^S,h = ( 
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(5) 



where we define x± = {Uo + Ep) sin P±/ko{hvp)^ in 
which fco is defined as ko = (( ^°+^^ )^ - (j2)i/2_ Eq. 



5|) e*'^* is defined as u±/v± and 
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In the case of d-wave symmetry, the orientational de- 
pendency of superconducting gap reads as A(7±) = 
A(r) cos(27± — 2a) in which a represents the orientation 
angle of d-wave superconducting gap. We now proceed 
and using the above Dirac-BdG wave functions and ap- 
propriate boundary conditions derive the Andreev and 



Normal reflection coefficients. By applying appropriate 
boundary conditions for the two interfaces which are lo- 
cated at a; = and L, we obtain all reflection and trans- 
mission coefficients. At last we assume that a large gate 
voltage Vb ^ 1 is applied into the narrow (i ^ 1) nor- 
mal region. In this case — VqL/IIvf is a constant which 
is called strength of barrier. The approximation is called 
thin barrier approximation regime in which the Normal 
region acts as an insulator. The Andreev and Normal 
reflection coefficients in the thin approximation regime 
are derived which are available in Appendix [A| for F|I|s- 
wave superconductor junctions. We assume a right-going 
electronlike quasiparticle within the ferromagnetic region 
incident into interface between the ferromagnetic and In- 
sulator regions, so the appropriate boundary condition in 
the interfaces at a; = is: 

and other boundary condition in interface between the 
insulator and superconductor regions ai x = L is: 

where rA,a and rjv,cr a-i'e amplitudes of spin-dependent 
Andreev and Normal reflection coefficients within the fer- 
romagnetic region, respectively. Other coefficients are 
transmission coefficients in the Normal and supercon- 
ducting regions. By applying the thin barrier approxi- 
mation on the obtained reflection and transmission fac- 
tors they reduce to simple factors that are available in 
the Appendix]^ For investigating the electronic thermal 
conductance of the junction one needs to calculate the 
probabilities of Andreev and Normal reflections namely 
Ij'A.g-P and |r7v,crp. By assuming a temperature gradient 
through the junction, the norm alized thermal conduc- 
tance r = r'/Fo is given as follo-wISIl^ 

f'/Fo =Y. / d€de„ cos{e„){l- I rMA^^^<y) ? 



(12) 



where Fq"^ = 2'K'^h?VFkB^a/ Ep is a constant. We pro- 
ceed to investigate the characteristics of electronic heat 
transport F of the mentioned junctions and throughout 
the paper we normalize energies with respect to Aq and 
we set Aq = fi = fc^ = 1 throughout our computations. 



III. ELECTRONIC THERMAL CONDUCTANCE OF THE 
F I S JUNCTIONS IN THE THIN BARRIER 
APPROXIMATION 

In this section we study electronic thermal transport 
characteristics of the Ferromagnetic! Insulator |s- wave su- 
perconductor junctions in the ballistic and thin barrier 




FIG. 2. (Color online) The normalized thermal conductance 
of F|I|S graphene-based junctions vs. normalized strength of 
barrier region Sl/vr for three values of h — 2Ao, 7Ao, 12Ao. 
The temperature has been fixed at T = 0.2Tc. 



approximation regime. In Fig. |2] we set T — 0.2Tc, 
and plot normalized thermal conductance F vs normal- 
ized strength of barrier JI/tt for three different values of 
magnetization texture strength h/Ao , also in Fig. [s] 
51 = is set and the normalized conductance is plot- 
ted for three values of temperatures vs magnetization 
texture strength /i/Aq. Throughout our calculations we 
have set Ep — lOAp and also used a large mismatch po- 
tential Uq. The normalized thermal conductance shows 
an oscillatory behavior vs f2/7r which this finding can be 
understood by noting the fact that how the amplitude of 
Andreev and Normal reflections depend on fi, (See Ap- 
pendix |A). In the thin barrier approximation, the width 
of normal layer L and barrier potential Vq set for small 
and large values, respectively. Andreev and Normal co- 
efficients are involved cos 257 and sin 2f2 terms which are 
periodic functions of Q and consequently the appeared 
periodic oscillations in the thermal conductance are orig- 
inated from the two appeared periodic functions in the 
Andreev and Normal coefficients. As it can be seen in 
Fig. [l] since the configuration contains one semi-infinite 
superconductor the Andreev bound states don't contri- 
bution to the tran sport ch aracteristics of the junction 
under consideratiorPEISlMl, Por small values of magnetic 
strength, the amplitude of oscillations has been enhanced 
in comparison with larger values of H/Aq. 

By increasing the magnetic strength, incident angle de- 
fined in Eq. (|3| reduce and hence diminish the amplitude 
of oscillations which means suppression of available prop- 
agating channels in the system. The normalized thermal 
conductance of the F|I|S junction is plotted vs. h/Ao, 
the magnetization strength of Ferromagnetic region, in 
Fig. |3]for three values of T = 0.2Tc, 0.5Ta, O.TTc and also 
r2 = is set for the three plots. The thermal conductance 
shows a minimum at /i ~ Ep and by increasing tempera- 
ture move the minimum towards smaller values of h. The 
magnetization texture splits Fermi level into two parts in 
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FIG. 3. (Color online) The normalized heat conductance 
of F|I|S graphene-based junctions vs magnetic exchange field 
strength h/Ao of Ferromagnetic region for three values of tem- 
peratures, T = 0.2Tc,0.5Tc,0.7Tc and fixed barrier strength 

at n = o. 



the fc-space and by increasing /i, the two parts separate 
upward and downward more and more. Increasing the 
exchange splitting suppress propagating Dirac Fermions 
modes in the configuration under consideration up to val- 
ues near h ~ Ep, for larger values of h the propagating 
channels enhance , see Ref.s [2] and UHl In the light of 
above discussion, the thermal conductance reach to its 
minimum value at /i ~ Ep that depends on the tempera- 
ture. The fact also can be inferred from Fig. [2] which the 
curve of h — 12Aq has an intermediate value between the 
curves oi h — 2Ao and h — TAq. We proceed to investi- 
gate effects of d-wave symmetry on the heat conductance 
of F|I|d junctions in the clean limit. 



IV. ELECTRONIC THERMAL CONDUCTANCE OF THE 
F]! D JUNCTIONS IN THE THIN BARRIER 
APPROXIMATION 

Now we present main results of this paper namely 
the fingerprints of d^2_y2-waye superconducting region 
on the electronic heat transport characteristics of F|I|d 
junctions in the thin barrier regime whose interfaces are 
located at x = 0, L. As it is seen in Fig. [T] unlike 
s-wave superconductors, the role of crystal orientation 
of (ij.2_j^2-wave superconductors with respect to inter- 
face is very important. We assume a two-dimensional 
d-wave superconductor with cylindrical Fermi surface 
in the /c-space is deposited on top of graphene sheet 
and connected to a sandwiched insulator region be- 
tween Ferromagnetic and superconducting regions. The 
pair potential for s-wave superconductor is isotropic 

I.e. A(r)=Ao tanh ^ 1. 76 a/Tc/T - 1. On the other 
hand, the pair potential for da;2_j,2-wave symmetry is 



T/T 

FIG. 4. (Color online) The normalized thermal conduc- 
tance r of F|I|d graphene-based junctions vs temperature 
for five values of d-wave superconducting gap orientation 
a = 0, 7r/16, 7r/8, 37r/16, 7r/4. The strength of magnetic ex- 
change field and barrier fixed at /i = 2Ao, $1 = 0, respectively. 



0-dependent, namely angle between the a- axis of the 
superconductor crystal and wavevector of the conduct- 
ing quasiparticles. In this case, superconducting gap is 
anisotropic i.e. A±{T,j)=Ad{T) cos{2j zL 2a) in which 
a is angle of a-axis with respect to normal trajectory 
to the interface (See Fig. [T]) and 7 is propagation an- 
gle of quasiparticles. The temperature dependency of 
d-wave superconductors is different from s-wave cas^^H, 
As it mentioned above, for the thin barrier approxima- 
tion is assumed that L <^ 1 and Vq 1, so one can 
consider 57 as a constant and terms involving fl reduce 
to simpler ones. Here we have set Ep = IOAq and use 
large mismatch potential Uq. Fig. |4] indicates electronic 
heat conductance of the F|I|d junctions vs temperature 
for five different values of crystal orientation of (1^2 _y2- 
wave superconductor a, the exchange field and strength 
of thin barrier are set a,t h ~ 2Ao and = 0, respec- 
tively. Electronic thermal conductance for a — 0,7r/16 
shows an exponential increase vs temperature and for 
larger values of a, the exponential form is modified to 
linear increase. The heat conductance shows completely 
linear increase vs temperature at maximum value of su- 
perconductor crystal orientation angle a — n/A that in- 
duces in mind the Wiedemann- Franz law for metals in 
low temperatures which thermal conductance is propor- 
tional to temperature, F cx T. The finding is arisen 
from orientational-dependent superconducting gap that 
increasing a decreases the propagating channels of super- 
conducting correlations described by Andreev refiection 
coefficients. Although the Dirac and Schrodinger equa- 
tions are used in graphene-based and metallic junctions 
respectively but the behaviors of thermal conductance 
in the graphene-based junctions are qualitatively similar 
to results of metallic N|I|d-wave junctions in which the 
propagating channels of moving quasi-particles is closed 
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FIG. 5. (Color online) The normalized thermal conductance 
r of F|I|d graphene-based junctions vs normalized strength of 
barrier region Q for seven values of superconducting gap ori- 
entation Q = 0,7r/16,7r/8,37r/16,3.37r/16,3.77r/16,7r/4. The 
strength of exchange field and temperature are fixed a,t h = 
2Ao, T = 0.2Tc, respectively. 
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FIG. 6. (Color online) The normalized thermal conductance 
of F|I|d graphene-based junctions vs strength of magnetic ex- 
change field h/ Aq for five values of superconducting gap ori- 
entation a = 0, 7r/16, 7r/8, 37r/16, 7r/4. The temperature and 
strength of barrier are fixed at T = 0.2Tc, $1 = 0, respectively. 



by increasing crystal orientation angle from to 7r/4 
in Ref.s and [5^1 In this context d-wave symmetry 
shows the same effects on thermal conductance of both 
graphene-based and metalhc junctions. The behaviors 
of heat conductance vs strength of barrier are shown in 
Fig. [5] for several values of a. Temperature and exchange 
field have set in h = 2Ao and T = 0.2Tc. The thermal 
conductance vs strength of barrier region shows an os- 
cillatory behavior and the increase of a enhance whole 
values of heat conductance. The period of oscillations vs 
strength of barrier suppresses completely for maximum 
crystal orientation angle a = n/i. In Fig. [6] the thermal 
conductance is plotted vs the strength of magnetization 
exchange field h/ Aq for several values of a and T — 0.2Tc. 
In general, T for F|I|d junctions vs h/Ao behaves similar 
to F|I|S configuration . The behavior can be verified by 
noting the mentioned reasons in the Sec. Ill for F|I|S 
case. Increasing the crystal orientation angle of d-wave 
superconductor up to a = 7r/4 can only enhance whole 
values of the thermal conductance vs h/ Aq. 



SUMMARY 



investigated as well. We found that for F|I|S junctions, 
the heat conductance vs magnetization strength H/Aq 
shows a minimum at values near h :^ Ep that by 
increasing temperature the minimum move towards 
smaller values of H/Aq. The finding is qualitatively 
similar to F|I|d junctions but increasing superconductive 
gap orientation a shifts whole values of T towards larger 
values and no change induces to the trend of F vs h/ Aq. 
The electronic thermal conductance vs barrier strength 
oscillates and shows identical behavior for F|I|S and F|I|d 
configurations for all values of superconductor crystal 
orientation a except values near a ~ 7r/4. By approach- 
ing to a = 7r/4, the propagating channels diminish and 
hence the amplitude of oscillations suppress. We found 
a Wiedemann- Franz law-like in the low temperature 
regime for thermal conductance of F\I\S junctions, 
namely F cx T. The electronic heat conductance shows 
an exponential growth vs temperature for small values 
of gap orientation angle a < 7r/8 and for larger values 
especially at a = 7r/4 approaches to completely linear 
growth, namely F cx T. 



In summary we have considered 
Ferromagnetic | Barrier | s /d-wstve superconductors 
graphene-Based junctions in the thin barrier approxima- 
tion and ballistic limit. We have utilized the Dirac-BdG 
equation and by employing Dirac-BdG wavefunctions 
derived the Andreev and Normal reflection amplitudes. 
Electronic thermal conductance F, of the two mentioned 
junctions in the thin barrier approximation has been 
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Appendix A: Andreev and Normal reflection coefFicients 
the thin barrier approximation regime for F|l S junctions 



Using the boundary condition Eq.s (10 11) and ap- 
plying the thin barrier approximation, the Andreev and 
Normal reflection coefficients for the F|I|S junctions are 
obtained as follows; 



rA 



rN 



y/cOS 6a \/cOS 6'^ COS 76 



Ti 



El = COS 2VL COS (- 



I sin P sin 7 — sin ( - 



-^)sin/3 



sin( 



cos /3 cos 7 + sin2f7cos(- 



Ti — cos (- 
T2 = cos (- 



I cos /3 cos 7 + sin 20. sin ( - 

Q ^ 

I sin j3 + cos 20, sin ( — 



-) sin P sin 7. 



The obtained coefficients recover the results of Ref.s S] 
[m and m for N|S, F|S and N|I|S graphene-based con- 
figurations, respectively. This can be justified by letting 
/i — >■ and — >■ in the above coefficients for F|I|S 
graphene-based junctions. 
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